Multi-Platform  Target  Detection  using 
Multi-Channel  Coherence  Analysis  and  Robustness 

to  the  Effects  of  Disparity 

Nick  Klausner*,  J.  Derek  Tucker*^,  and  Mahmood  R.  Azimi-Sadjadi* 

*  Department  of  Electrical  and  Computer  Engineering 
Colorado  State  University,  Fort  Collins,  Colorado  80523-1373 
Email:  {nklausne,dtucker,azimi}  @  engr.colostate.edu 
'  Naval  Surface  Warfare  Center  -  Panama  City  Division 
Panama  City,  FL  32407-7001 
Email:  james.d.tucker@navy.mil 


Abstract — The  use  of  multiple  disparate  platforms  in  many 
remote  sensing  and  surveillance  applications  allows  one  to  exploit 
the  coherent  information  shared  among  all  sensory  systems 
thereby  potentially  reducing  the  risk  of  making  single-sensory  bi¬ 
ased  detection  and  classification  decisions.  This  paper  introduces 
a  target  detection  method  based  upon  multi-channel  coherence 
analysis  (MCA)  framework  which  optimally  decomposes  the 
multi-channel  data  to  analyze  their  linear  dependence  or  coher¬ 
ence.  This  decomposition  then  allows  one  to  extract  MCA  features 
that  can  be  used  to  implement  a  coherence-based  detector.  This 
detector  is  applied  to  a  data  set  of  simulated  disparate  sonar 
imagery  provided  by  the  Naval  Surface  Warfare  Center  (NSWC) 
-  Panama  City.  This  database  contains  images  of  both  targets  and 
non-targets  with  various  variabilities  with  respect  to  resolution, 
signal-to-noise  ratio  (SNR),  target  and  non-target  types,  etc. 
Sensitivity  analyses  are  then  carried  out  in  order  to  gauge  the 
performance  under  such  variablities  that  may  be  encountered 
in  disparate  multi-platform  detection  problems.  Performance  of 
the  detection  method  will  be  given  in  terms  of  probability  of 
detection  (Pd),  probability  of  false  alarm  (P/a),  and  the  receiver 
operating  characteristic  (ROC)  curves. 

Index  Terms — Binary  hypothesis  testing,  disparate  sonar  plat¬ 
forms,  multi-channel  coherence  analysis,  underwater  target  de¬ 
tection 

I.  Introduction 

The  problem  of  developing  a  robust  underwater  target  de¬ 
tection  and  classification  system  that  can  operate  with  multiple 
disparate  sensing  systems  poses  many  technical  challenges 
and  has  recently  attracted  considerable  interest  in  many  areas 
of  remote  sensing.  In  the  traditional  centralized  processing 
scheme  preliminary  detection,  feature  extraction,  and  object 
classification  are  performed  based  upon  the  data  collected 
by  every  individual  sensor  platform.  Final  decision  making 
usually  takes  place  at  the  fusion  center  using  some  type  of 
decision-level,  feature-level,  or  combined  fusion  mechanism. 
However,  decision-making  based  upon  individual  sensory  data 
can  lead  to  incomplete,  degraded,  or  locally  biased  (sensor- 
level)  decisions  resulting  in  a  poor  final  detection  and  classi¬ 
fication  performance  at  the  fusion  center.  Thus,  new  methods 
are  needed  to  detect,  isolate,  and  represent,  in  terms  of 


some  pertinent  attributes,  the  coherent  or  common  information 
shared  among  the  multiple  data  sets.  This  can  be  an  extremely 
challenging  problem  due  to  the  disparate  nature  of  the  problem 
and  variations  in  the  operating  conditions  which  can  often  arise 
in  many  multi-platform  sensory  systems. 

An  approach  to  multiple-channel  signal  detection  was  con¬ 
sidered  in  [1]  by  introducing  the  Generalized  Coherence  (GC) 
estimate,  a  nonparametric  measure  of  the  information  shared 
by  any  number  of  channels.  This  measure  is  formed  by  finding 
the  determinant  of  the  Gram  matrix  formed  by  all  N  channels 
normalized  by  the  product  of  the  squared  norm  of  each 
channel.  This  GC  measure  exhibits  many  attractive  properties 
including  invariance  to  constant  channel  gain.  The  GC  measure 
is  shown  to  have  an  appealing  geometrical  interpretation  as  it 
is  noted  that  the  determinant  of  the  Gram  matrix  heuristically 
measures  the  squared  volume  of  a  parallelepiped  in  n-D 
space,  where  n  corresponds  to  the  dimension  of  each  channel. 
Thus,  the  more  the  channels  have  in  common,  the  smaller  the 
normalized  volume,  and  the  closer  the  GC  estimate  to  one. 
The  distribution  of  the  GC  estimate  is  then  derived  under 
conditions  of  white  Gaussian  noise,  theoretically  allowing 
one  to  determine  thresholds  corresponding  to  a  particular 
false  alarm  rate.  Simulations  are  then  provided  demonstrating 
the  GC  estimate’s  capability  in  multiple-channel  detection 
problems.  The  existing  work  [2]  -  [3]  in  the  area  of  target 
detection  from  sonar  imagery  has  primarily  been  focused  on 
one  sonar  platform,  with  fusion  across  multiple  algorithms.  In 
[4],  the  adaptive  clutter  filter  detector  in  [3]  is  individually 
applied  to  three  different  sonar  images  varying  in  frequency 
and  bandwidth.  Final  classification  is  done  using  an  optimal 
set  of  features  using  a  nonlinear  log-likelihood  ratio  test  where 
the  decisions  of  the  individual  detector  and  classifier  are 
fused.  The  optimal  set  of  features  is  determined  based  upon 
cascading  another  classifier  on  the  previous  classifier  during 
the  training  stage.  This  is  done  as  a  repeated  application  during 
the  training  stage  where  at  each  iteration  the  threshold  and 
optimal  feature  set  is  chosen  and  updated.  In  our  previous 
work  [5],  we  developed  a  new  framework  for  multi-sensor 
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coherence  analysis  using  Multi-Channel  Coherence  Analysis 
(MCA)  [6]  that  can  be  applied  to  the  data  collected  using 
multiple  disparate  sonar  systems.  Using  this  method  allows  for 
the  simultaneous  detection  and  feature  extraction  of  coherent 
target  information  among  multiple  disparate  sonar  images. 

This  paper  reviews  the  A’ -channel  coherence-based  detector 
using  the  MCA  framework  [5].  This  detector  exploits  the 
coherence  of  objects  present  in  N  disparate  channels  based 
on  the  assumption  that  the  presence  of  objects  in  all  data  sets 
will  lead  to  a  higher  level  of  coherence  compared  to  that  of 
noise  alone.  New  expressions  for  the  log-likelihood  ratio  and 
J-divergence  in  the  MCA  framework  are  provided  and  used  for 
the  simultaneous  detection  of  targets  from  N  disparate  sonar 
data.  In  order  to  examine  the  proposed  method’s  robustness 
to  variables  of  disparity,  a  thorough  sensitivity  analysis  was 
conducted  using  a  data  set  consisting  of  synthetically  gener¬ 
ated  images  of  targets  provided  by  the  NSWC  -  Panama  City 
Division. 

The  characterizing  difference  among  the  method  presented 
in  this  work  and  that  of  [1]  are  the  statistics  used  to  form 
detection  decisions  on  whether  or  not  an  unknown  but  common 
signal  is  present  within  the  data  collected  by  a  number  of 
sensing  systems.  While  the  method  presented  in  [1]  forms  the 
inner  product  of  each  pair-wise  combination  of  channels  to 
measure  the  overall  coherence,  the  method  of  this  work  first 
projects  the  multi-channel  data  into  a  new  coordinate  system 
where  the  relations  (linear)  among  the  channels  are  maxi¬ 
mized.  The  optimal  Gauss-Gauss  detector  [7]-  [8]  for  binary 
hypothesis  testing  in  the  composite  multi-channel  domain  is 
then  developed  by  deriving  expressions  for  the  log-likelihood 
ratio  and  J-divergence,  a  global  measure  of  the  amount  of 
discriminatory  information  available  for  detection.  However, 
one  advantage  of  the  statistic  presented  in  [1]  is  the  fact  that 
the  measure  is  nonparametric  in  the  sense  that  no  a  priori 
information  of  the  problem  is  needed  before-hand. 

This  paper  is  organized  as  follows:  Section  II  will  review 
the  MCA  framework.  Section  III  develops  the  MCA-based 
Gauss-Gauss  detection  method  [7]-  [8].  Section  IV  provides 
the  results  and  discussion  of  the  sensitivity  analyses  and  finally 
concluding  remarks  will  be  made  in  Section  V. 

II.  Multi-Channel  Coherence  Analysis 

Consider  N  zero  mean  random  vectors,  xi,  X2,...,  and  xy, 
representing  multiple  data  channels  comprising  the  composite 
data  channel  z  =  [x^  x^  •  •  •  xj^]H  £  Cdxl.  Without  loss 
of  generality,  we  will  assume  all  random  vectors  to  be  zero 
mean  throughout  this  analysis.  Let  each  channel  Xj  £  Cdj  Xl 
be  of  dimension  dj,  where  it  is  assumed  that  xi  is  of  the 
smallest  dimension  and  we  denote  d  =  dj-  The  d  x  d 

dimensional  covariance  matrix  of  the  composite  data  channel 
z  is  given  by 

-Rn  R12  ■  ■  ■  Rin 

R21  R22  ■  ■  ■  R2N 

.  .  ._  .  ,  (1) 

Rni  Rn  2  ■  ■  ■  RnN 


where  Rjk  =  E[xjxf?]  is  the  auto-covariance  (j  =  k )  or 
cross-covariance  (j  ^  k )  matrices  of  data  channels  Xj  and  x^ 
and  clearly  we  have  Rjk  =  Rkj- 

Similar  to  Canonical  Correlation  Analysis  (CCA)  [9],  [10] 
the  ith  multi-channel  coordinate  of  the  jth  channel  is  found 
by  searching  for  the  ith  coordinate  mapping  vector,  ayj,  of 
data  channel  Xj.  This  linear  transformation  produces  the  ith 
multi-channel  coordinate  for  the  jth  channel, 

Vij  =  Oi^jXj.  (2) 

If  the  ith  coordinate  mapping  vectors  are  found  for  all 
N  channels,  we  can  then  obtain  the  composite  coordinate 
mapping  vector  a,  =  af2  ■  ■  ■  afN]  which  is 

then  used  to  find  the  composite  coordinate  vector  v,  = 

Ki  W,2  •  •  •  Wn]T  =  [xf  a^i  xf  •  •  •  x^aitN] 
which  consists  of  the  ith  multi-channel  coordinate  of  every 
channel.  The  associated  covariance  matrix  of  vt  is  given  by 

RviVi  =  E  [vjvf  ] 

ai^lRllai,l  '  '  '  OlfpRlNOli'N 

:  ••.  ;  (3) 

■  '  '  OcfNRNNOLi^N 

Recall  that  in  the  two-channel  CCA  [7],  [11],  the  correla¬ 
tions  between  the  mapped  coordinates  are  maximized  subject 
to  the  constraint  that  the  transformed  coordinates  have  unit 
variance.  In  the  multi-channel  case,  however,  the  analysis  is 
not  as  well-defined  as  all  correlations  between  all  possible 
pairs  of  channels  must  be  maximized  simultaneously.  To 
accomplish  this,  one  approach  that  has  been  proposed  [6]  is  to 
maximize  the  sum  of  all  correlations  (the  SUMCOR  objective 
function)  subject  to  the  unit  trace  constraint  of  matrix  RViVi. 
Thus,  the  optimization  problem  for  finding  the  first  composite 
coordinate  mapping  vector  a  t  using  the  objective  function  and 
constraint  just  described  becomes 

N  N 

ai  =  arg  max  ^  ^  a^dRjtkaltk 
j= 1  k= 1 
N  N 

=  arg  max  EE  (4) 

j= 1 fc= 1 

subject  to  the  constraint 

N 

tr (-RviVl)  =  'y  ,ai,jRjjai,j  =  1j 

j- 1 

It  is  shown  [6]  that  the  constrained  optimization  problem  for 
the  first  coordinate  mapping  vectors,  aij  using  a  Lagrange 
multiplier  method  leads  to 

N 

y  'j  kij cr  1  i;  —  XiRjjaij,  V  j,  k  £  [l,iV] 

fc= 1 

or  in  matrix  notation  as 


-7?zz^i  —  Ai-Z7a.r, 


(5) 


where  D  is  a  block  diagonal  matrix  with  diagonal  blocks 
Rjj,  V  j  G  [1,7V],  i.e. 

D  =  diag  [Ru,R,22,  ■  ■  ■  ,Rnn]  •  (6) 

The  result  in  (5)  represents  a  generalized  eigenvalue  problem 
for  which  standard  methods  of  solution  are  well-known.  We 
will  then  consider  the  simultaneous  solution  to  all  mapping 
vectors  a^s,  i  G  [l,d]  and  write  (5)  as  RZZA  =  DAA 
where  A  consists  of  all  d  coordinate  mapping  vectors,  and 
A  consists  of  all  d  eigenvalues.  This  solution  can  then  be 
rewritten  in  terms  of  a  standard  eigenvalue  problem  EP  = 
PA  where  E  =  D~^RZZD~^  and  P  is  a  unitary  matrix 
(ppH  =  pH p  _  /  j  clearly,  we  can  find  the  mapping  matrix 
A  via  A  =  D~^P.  Inspection  of  matrix  E  shows  that  it 
is  simply  the  composite  covariance  matrix  of  the  whitened 
version  of  z  =  [  ]H .  That  is,  if  we  define 

this  whitened  version  of  the  composite  data  channel  vector 
by  w  =  [  wf*  •  •  •  ]H  =  D~ 2 z,  then  the  whitened 

composite  vector  w  has  correlation  matrix  E  [ww;/]  = 

D~^  RZZD~^  =  E.  Matrix  P  is  then  used  to  map  the 
whitened  channels  to  their  multi-channel  coordinates.  In  or¬ 
der  to  find  mapping  vectors  corresponding  to  the  principal 
coordinates  [6],  we  only  consider  the  r  =  di  =  min;  {dj} 
coordinates  such  that  Ai  >  A2  >  ...  >  Ar.  Thus,  A  = 
diag  [Ai,  A2, . . . ,  Ar]  will  become  a  r  x  r  diagonal  matrix 
composed  of  the  dominant  eigenvalues  and  P  will  become  a 
dxr  matrix  composed  of  the  eigenvectors  corresponding  to  r 
dominant  eigenvalues.  To  find  the  mapped  coordinate  vector, 
v,  that  contains  all  mapped  coordinates  for  all  TV  channels, 
we  will  first  define  T ' j  (dimension  dj  x  r)  to  contain  those 
dominant  r  eigenvectors  Pij,  V  t  G  [l,r]  of  the  mapping 
matrix  P  that  correspond  to  the  jth  channel 

=  [  Pl,j  P2 ,j  ■  ■  ■  Pr,j  ]  ,  V  j  G  [1,  TV]  .  (7) 

Clearly,  the  connection  between  P  and  is  evident 

^2 

P  =  .  (8) 

-  J  dxr 

Note  that  in  the  case  of  two  channels,  'k  1  and  T,  are  directly 
related  to  the  mapping  matrices  of  CCA  [8].  All  of  the  mapped 
coordinates  of  the  jth  channel  can  then  be  found  by 

//.,  -  Xr  V  je  [1,7V],  (9) 

where  fij  =  [  Vij  V2j  ■  ■  ■  vrj  }T .  Clearly,  we  have  the 
following  two  properties 

N  N 

3= 1  3= 1 

N  N  N  N 

EE£W1  =  EE$f^w^  =  a 

j= 1  k—  1  j= 1  k—1 


If  we  define  block  diagonal  matrix  Sk  that  contains 
the  T'j  matrices  along  its  diagonal  blocks,  i.e.  fir  = 
diag  [^7r,  ’f'2,  •  ■  • ,  ^jv],  then  we  can  resolve  all  TV  channels 
into  their  multi-channel  coordinates  using 

v  =  A/hw  =  A>h  (10) 

III.  MCA  Detection 

A  classical  detection  problem  is  that  of  choosing  between 
two  hypotheses  that  are  relevant  to  the  given  problem.  For 
this  coherence-based  detector,  the  null  hypothesis  (TTo)  is  the 
hypothesis  that  all  TV  channels  consist  of  background  noise  and 
the  alternative  hypothesis  (Pi)  that  all  TV  channels  consist  of 
signal  plus  noise.  Figure  1  shows  the  graphical  setup  of  the 
problem  under  consideration.  Several  simplifying  but  sensible 
assumptions  used  in  this  analysis  are 

1)  Noise  between  different  channels  is  mutually  uncorre¬ 
lated,  i.e.  E  [rijn^]  =  0  V  j,  fc  G  [1,  TV],  j  ^  k. 

2)  Signal  is  uncorrelated  with  the  background  noise,  i.e. 
E  isink  ]  =  E  [nisf]  =  0  V  j,  k  G  [1,  TV], 

3)  Noise  contained  on  any  one  channel  has  covariance 
matrix,  i.e.  E  [11,-11^]  =  Rn.  V  j  G  [1,7V], 

4)  Signal  contained  on  any  pair  of  channels  has  covariance 
matrix,  i.e.  E  [sj-sf  ]  =  RSjk  V  j,  k  G  [1,  TV]. 


Fig.  1.  Multi-Channel  Hypothesis  Test  and  MCA. 


Under  //q,  the  matrices  Rzz  and  D  become 

Pzz0  E0  —  diag  [Pni ;  T?n2 ,  •  •  •  ,  PriAr]  ■ 

Note  that  the  subscript  notation  refers  to  the  hypothesis  being 
considered.  Thus,  the  solution  to  the  eigenvalue  problem  under 
Po  leads  to  Ao  =  /,  while  for  .4q  no  unique  solution  exists. 
Now,  since  any  non-zero  vector  will  satisfy  the  generalized 
eigenvalue  problem,  for  simplicity  we  choose  the  eigenvectors 
of  the  null  hypothesis  to  be  the  same  as  those  of  the  alternative 
hypothesis,  i.e.  Ao  =  A\  =  A. 

Under  Pi  and  using  the  stated  assumptions,  the  correspond- 


ing  Rzz  and  D  matrices  are 

RS11  d-  ^ni  Rsi2 

Rsin 

Rs21  Rs22 

R&2N 

Rzzi  —  .  . 

RsN1  R&N2 

Rsnn  A  RnN 

Di  diag  [E£Sll  H-  Rni ,  Rs22  Rn 2  5  • 

■  ■  J  Rsnn  d"  R^n  ] 

leading  to  the  following  arbitrary  eigenvalue  problem. 
Pzzi  Ai  =  D1A1A1 


(ID 


The  log-likelihood  ratio  that  minimizes  the  risk  involved  in 
deciding  between  the  two  hypotheses  is  defined  [12]  to  be 


Rearranging  and  using  the  cyclic  property  of  the  trace,  we  can 
write  the  J-divergence  as 


l(z)  =  In 


Ppl-Hi) 

p(z\Ho) 


(12) 


which  is  compared  to  a  global  threshold  rj  to  form  the 
discriminant  function  'y(z). 


7(z)  =  /  1  ~  Hl  1 (z)  ^  ?? 

71 ’  \  0  ~  H0  l{ z)  <  V 

Assuming  that  under  both  hypotheses  the  composite  data  chan¬ 
nel  z  is  multivariate  Gaussian  with  zero  mean  and  covariance 
matrix  i?zz,  the  log-likelihood  ratio  of  the  composite  data 
vector  becomes 


l{ z)  =  (Rj0  -  RZZl)  z  (13) 

Next,  we  can  formulate  R~z  in  terms  of  the  sum  of  the 
correlations  of  each  coordinate  and  their  corresponding  eigen¬ 
vectors.  To  do  this,  we  recall  the  fact  that  PH EP  =  A. 
Taking  the  inverse  of  this  relationship,  it  is  simple  to  show 
that  Rfzl  =  AA~1AH .  Thus,  the  log-likelihood  function  in 
(13)  becomes 

l{ z)  =  zH  (A0Aq1Aq  -  AtA^A? )  z. 

Since  the  eigenvalues  under  the  null  hypothesis  are  all  one, 
i.e.  Ao  =  /,  and  owing  to  the  lack  of  a  unique  solution  for 
the  mapping  matrices  of  the  null  hypothesis,  we  can  write  the 
log-likelihood  ratio  as 

Kz)  =  zH  [Ai  (I  -  A]-1)  Af  ]  z,  (14) 

where  A\  and  Ai  are  the  mapping  matrix  and  diagonal  matrix 
of  multi-channel  correlations,  respectively,  for  the  set  of  data 
with  which  we  are  performing  the  hypothesis  test. 

Next,  we  will  formulate  the  J-divergence  [11],  [13]  which 
is  a  measure  of  the  separability  of  the  two  hypotheses.  The 
J-divergence  is  defined  to  be 

J  =  EHl[l(z)\-EHo[l(z)\,  (15) 

where  Eh1  [  •  ]  and  EH[:  [  ■  ]  represent  the  expectation  opera¬ 
tion  evaluated  under  the  Hi  and  Hq  hypotheses,  respectively. 
The  expected  value  of  the  log-likelihood  function  becomes 

E  [l(z)\  =  E  [tr  (zH Qz)]  .  (16) 

Where  Q  =  (  RfJ()  —  RZZl )  •  Using  the  cyclic  property  of  the 
trace,  we  can  write 

E[l(z)\  =  E  [tr  (Qzzh)\ 

=  tr  (QRZZ)  ■  (17) 

Thus,  we  can  write  the  J-divergence  as 

J  =  EHl  [Z( z)]  -  EHo  [Z( z)] 

=  tr  (QRZZl )  -  tr  ( QRZZo ) 

=  tr  [—2/  +  Rj0RZZl  +  RZZ1RZZ0\  ■  (18) 


J  —  tr  [ — 2 1  +  Ai  +  A-l  —  ''y  ^  (— 2  +  A,;  +  A.j  .  (19) 

i— 1 

Therefore,  the  only  pieces  of  information  we  need  to  know 
when  performing  detection  in  this  framework  are  the  matrices 
A\  and  Ai. 


IV.  Simulation  Results 

A.  Data  Description  and  Pre-Processing 

This  (V-channel  coherence-based  detection  method  was  then 
applied  to  a  data  set  consisting  of  synthetically  generated  sonar 
images  (snippets)  of  both  targets  and  non-targets  of  different 
geometrical  shapes.  The  sonar  snippets  were  generated  with 
different  resolutions,  SNR,  range,  and  aspect  angles  mimicking 
different  realistic  operating  conditions.  For  this  study,  two 
different  resolutions,  namely  1  vn  and  3 in,  were  considered. 
Additionally,  SNR  ranged  from  0 dB  to  15 dB  in  increments  of 
3 dB,  range  values  spanned  from  10m  to  120?n  in  increments 
of  lrn,  and  aspect  angle  ranged  anywhere  from  0°  to  360° 
in  increments  of  1°.  A  subset  of  the  data  corresponding 
to  target  snippets  was  used  to  represent  the  Hi  hypothesis 
while  background  snippets  were  used  to  represent  the  Hq 
hypothesis.  Thus,  all  non-target  snippets  were  excluded.  The 
subset  of  target  snippets  (1610  snippets)  was  further  parti¬ 
tioned  into  3  different  parts  forming  138  cone-type  targets,  736 
cylinder-type  targets,  and  736  trapezoidal-type  targets.  When 
performing  detection,  each  of  the  N  images  is  partitioned 
into  blocks  of  size  dependent  on  the  resolution.  Each  block 
is  then  channelized  and  used  to  form  a  realization  of  the 
composite  data  vector,  z.  The  ensemble  set  formed  from  all 
realizations  pertaining  to  that  set  of  N  images  is  then  used  to 
form  the  sample  composite  covariance  matrix  /i'zz  which  is 
then  decomposed  by  the  MCA  process  and  used  to  form  the 
discriminant  function,  l(z)  in  (14).  Each  block  (realization) 
from  the  ensemble  set  is  then  applied  to  the  log-likelihood 
ratio  test  and  if  50%  or  more  of  the  blocks  pass,  that  set  of 
N  snippet  images  is  said  to  contain  a  target.  Figure  2  displays 
the  process  behind  the  MCA-based  detection  system  used  for 
this  analysis. 

B.  Dual  Resolution  Disparate  Detection  Performance 

In  some  disparate  detection  applications,  each  platform 
may  carry  multiple  sensing  systems  with  different  spatial  and 
spectral  characteristics  in  order  to  highlight  different  attributes 
of  the  target.  To  simulate  such  a  situation,  a  two-channel  de¬ 
tection  problem  was  constructed  where  each  channel  consisted 
of  snippets  of  targets  of  the  same  type  at  the  same  range, 
aspect  angle,  and  SNR.  However,  the  two  channels  differed  in 
resolution,  one  snippet-image  of  high  resolution  (1  in)  and  the 
other  of  a  lower  resolution  (3 in).  When  performing  detection, 
a  4  x  4  block  size  was  used  for  the  high  resolution  images  and 
a  2  x  1  block  size  for  those  of  lower  resolution.  This  setup 
was  then  run  for  all  1610  images  at  various  ranges  and  aspect 


Image  1 


Fig.  2.  N-Channel  MCA  Detector. 


angles  and  the  results  partitioned  on  the  basis  of  target  type 
and  SNR. 

Figures  3(a)-(c)  display  the  ROC  curves  for  conical,  cylin¬ 
drical,  and  trapezoidal  type  targets,  respectively,  at  three 
different  SNR  values  of  0 dB,  6dB ,  and  12 dB,  respectively. 
Figures  4(a)-(c),  on  the  other  hand,  display  sensitivity  to  target 
types  for  a  given  SNR.  Table  I  gives  the  probability  of  detec¬ 
tion  at  the  knee  point  of  the  ROC,  i.e.  where  +  Pfa  =  1, 
for  all  target  types  and  SNR  values.  As  can  be  observed 
from  the  results  in  Table  I,  it  is  apparent  that  for  cylindrical 
and  trapezoidal  targets  the  detection  performance  generally 
improves  as  a  function  of  SNR  as  one  would  expect.  However, 
the  performance  of  the  detector  for  the  cone  targets  does  not 
follow  the  same  behavior.  This  may  be  attributed,  in  part,  to 
the  fact  that  only  a  small  number  of  cone  targets  were  available 
for  this  study. 

TABLE  I 

Probability  of  Detection  (%)  vs.  SNR 


Target  Type 

0  dB 

3  dB 

6  dB 

9  dB 

12  dB 

15  dB 

Cone 

91.30 

94.93 

96.38 

89.86 

81.16 

87.68 

Cylinder 

83.70 

85.19 

82.20 

85.33 

89.67 

94.02 

Trapezoid 

84.24 

84.51 

84.65 

85.73 

90.08 

93.75 

C.  Dual  Aspect  Angle  Separation  Disparate  Detection  Perfor¬ 
mance 

Another  question  that  may  arise  in  any  disparate  detection 
problem  is  that  of  sensor  locations.  That  is,  for  different  types 
of  targets  how  does  the  detection  performance  change  as  a 
function  of  disparities  in  location  of  the  sensor  platforms.  This 
clearly  relates  to  the  target’s  aspect/orientation  with  respect  to 
the  sensor  platform  as  well  as  the  range  of  the  platform  from 
the  target.  To  determine  the  answer  to  this  question,  a  study 
is  carried  out  where  both  channels  consisted  of  images  of  the 
same  resolution  (1  in),  at  ranges  within  ±1to  of  one  another, 
and  at  an  identical  SNR  of  9 dB,  while  the  disparateness  was 
with  respect  to  aspect  angle  separation.  More  specifically, 
the  two  channels  correspond  to  sonar  snippet-images  from 
the  same  target  at  two  aspect  angles  with  separation  angle  9 
such  that  if  0i  and  0 2  represent  the  aspect  angles  associated 
with  their  respective  image  then  pairs  of  images  were  chosen 


such  that  |0i—02|  G  [0  —  5,9  +  8\,  where  6  represents  the 
perturbation  from  the  separation  angle,  9,  due  to  non-uniform 
motion  of  the  vehicle.  Here,  the  value  of  5  was  chosen  to  be 
10°.  The  aspect  angle  separation  9  was  then  varied  from  0°  to 
180°  in  increments  of  30°  and  its  affect  on  the  performance  of 
the  detector  was  studied.  Figures  5(a)  and  (b)  display  the  ROC 
curves  for  cylindrical  and  trapezoidal  targets,  respectively,  for 
several  values  of  9.  Note  that  all  images  in  the  database 
corresponding  to  cone  targets  were  generated  at  an  angle 
of  0°  and  thus  excluded  from  this  study.  Table  II  gives  the 
probability  of  detection  at  the  knee  point  of  the  ROC  versus 
aspect  angle  separation  and  target  type.  From  both  Figures  5(a) 
and  (b)  and  the  results  in  Table  II,  it  can  be  concluded  that  the 
performance  of  the  detector  is  fairly  robust  to  disparateness  in 
aspect  angle  separation  as  the  probability  of  detection  at  the 
knee  point  of  the  ROC  never  falls  below  92%. 


TABLE  II 

Probability  of  Detection  (%)  vs.  Aspect  Angle  Separation  (8) 


Target  Type 

0° 

30° 

O 

O 

CD 

90° 

120° 

150° 

180° 

Cylinder 

94.40 

96.14 

97.83 

98.33 

97.47 

96.82 

94.92 

Trapezoid 

92.86 

96.29 

97.21 

97.66 

97.57 

96.09 

93.14 

Figure  6  displays  the  plots  of  empirical  J-divergence  as  a 
function  of  aspect  angle  separation,  9,  for  both  cylindrical  and 
trapezoidal  targets.  Both  curves  were  generated  by  empirically 
estimating  the  difference  in  means  of  the  log-likelihood  ratio 
among  target  {Hi)  and  noise  {H0)  snippets  and  averaging 
over  all  such  pairs  of  images  that  match  the  criteria  explained 
previously  (i.e.  same  SNR,  same  resolution,  ranges  within 
lm,  aspect  separation  within  some  range  of  a  particular 
angle).  The  results  in  Figure  6  match  what  was  observed  in 
Figures  5(a)  and  (b)  as  the  performance  of  the  detector  seems 
to  improve  as  the  separation  in  aspect  angle  deviates  from 
0°  and  180°.  Again  we  can  draw  the  same  conclusion  that 
the  detector  is  fairly  robust  to  separation  in  aspect  angles  as 
the  difference  between  maximum  and  minimum  J-divergence 
never  grows  larger  than  approximately  0.6.  Additionally,  the 
J-divergence  values  for  cylindrical  target  are  higher  than  those 
of  trapezoidal  target.  This  result  is  also  evident  in  the  ROC 
plots  of  Figures  5(a)  and  (b). 


(a)  Cone  Target  Type 
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(b)  Cylindrical  Target  Type 


(b)  6dB  SNR 


(c)  Trapezoidal  Target  Type 


Fig.  3.  Detector  Performance  for  Different  Target  Types  vs.  SNR.  Fig-  4.  Detector  Performance  for  Various  SNR  Values  vs.  Target  Type. 


V.  Conclusion 

In  this  paper,  a  multi-channel,  multi-sensory  binary  hypoth¬ 
esis  detection  system  has  been  introduced  using  the  MCA 
framework.  An  A' -channel  Gauss-Gauss  detector  was  then 
formulated  in  the  MCA  coordinates.  Detection  is  performed 
by  extracting  the  multi-channel  mapping  vectors  and  the  cor¬ 


relation  sums  from  the  data  samples  collected  by  the  sensory 
platforms.  These  mapping  vectors  and  coherent  features  are 
then  used  in  the  log-likelihood  ratio  to  detect  targets  in  the 
sonar  images.  This  MCA-based  detector  is  then  applied  to  a 
synthetic  sonar  image  data  set  to  investigate  the  performance 
of  the  detector  to  different  sensing  disparities.  These  studies 
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(a)  Cylindrical  Target  Type 


(b)  Trapezoidal  Target  Type 


Fig.  5.  Detector  Performance  for  Two  Target  Types  vs.  Aspect  Angle 
Separation. 


were  specifically  focused  on  the  dual  channel  detection  prob¬ 
lem  as  the  data  set  did  not  allow  for  more  than  2  channels. 
We  evaluated  the  detector  performance  versus  differences 
in  target  types,  SNRs,  and  aspect  angles.  The  detector  was 
found  to  be  robust  to  SNR  changes  across  most  of  the  target 
types  and  robust  to  varying  degrees  of  aspect  separation  by 
observing  the  probability  of  detection  at  the  knee  point  of 
the  ROC  curve.  Furthermore,  the  J-divergence  was  empirically 
estimated  and  shown  to  also  be  robust  to  disparities  in  aspect 
angle  separation.  Through  this  work  we  have  shown  that  MCA 
provides  an  elegant  framework  when  performing  coherence- 
based  detection  among  multiple  disparate  sensory  channels. 
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